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Detection of Ceo in the proto-planetary nebula IRAS 01005+7910 
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ABSTRACT 

We report the first detection of buckminsterfullerene (Ceo) in a proto- 
planetary nebula (PPN). The vibrational transitions of Ceo at 7.0, 17.4, and 
18.9 /xm are detected in the Spitzer/IRS spectrum of IRAS 01005+7910. This 
detection suggests that fullerenes are formed shortly after the asymptotic giant 
branch but before the planetary nebulae stage. A comparison with the obser- 
vations of Ceo in other sources is made and the implication on circumstellar 
chemistry is discussed. 



Subject headings: infrared: stars 
stellar matter 



stars: AGB and post-AGB — stars: circum- 



1. Introduction 

Fullerenes, together with other forms of carbon such as graphite, dia monds, and car- 



byne s, are expected to be a important components of the interstellar medium (jHenning fc Salama 



19981 ). Th e most stable fulle rene is buckminsterfullerene (Cqq) which has a soccer-ball like 



interstellar bands (e.g. 



standing mystery (see iHerbig 



Kroto et al. 



1995 



Lesei et al. 


1988), the origin c 


review) . 


Foine & Ehrenfreund 



the laboratory spectrum of C^~ in argon and neon matrices shows approximate matches 
with two diffuse interstellar bands in the near infrared spectra of stars. Possible formation 
processes of fullerenes in space that have previously been discussed include condensation in 
supernova gas, shock-induced decomposi tion of hydrogenated amorphous carbon g rains, cold 



interstellar gas-phase chemistry, etc (see iMoutou et al.lll999l : ISellgren et al.ll2010l and refer- 
ences there in). The carbon-rich hydrogen-poor circumstellar envelopes, such as Wolf-Rayet 
(WR) stars and R Cor onae Borealis stars, have also been proposed as favorable sites for the 



synthesis of fullerenes (jCherchneff. et al.ll2000i : iGoeres et al.lll992l ). However, none of these 



theoretical predictions have been confirmed observationally. 
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The search for fullerenes in space started soon after their synthesis in the laboratory. 
The early efforts to search for the electronic transitions of Cgo in the optical and UV wave- 
lengths have resulted in no definite evid ence showing its presence in space (e.g. lSnow fc Seab 



19891 ; iHerbig Il2000l ; ISassara et al.ll200ll). Another a p proach is to s earch for the infrared vibra- 
tional modes of the molecule (jClayton et al.lll995l ). iKwok et al.l (119991 ) noted a faint feature 
at 17.85 and 18.90 //m in the ISO-SWS spectrum of the proto-planetary nebula (PPN) IRAS 
07134+1005 and suggested that they could correspond to the z/ 2 7 and i/ 28 vibrational modes 
of the C 60 molecule. However, these features are not confirmed by subsequent Spitzer obser- 
vations. A search of these ban ds in the reflection nebula NGC 7023 with ISO was also not 
successful (IMoutou et al.lll999). Only very rece ntly, C 60 was unambiguously detected in the 



PNs (IGarcia-Hernandez et al.ll2010l ) 



planetary n ebula (PN) Tc 1 (|Cami et al.ll2010f ) and the reflection nebulae NGC 7023 and 
NGC 2023 (pellgren et al.l 120101) . T his was followed by the detections of Cqq in four more 



These detections of Ceo in objects at the late stages of stellar evolution raise the con- 
crete question of how Ceo is formed in circumstellar environment. Since the element C is 
synthesized in Asymptotic Giant Branch (AGB) stars and many C-based molecules have 
been detected in the outflow of AGB stars, it is reasonable to expect the C 60 can also be 
one of the products of circumstellar molecular synthesis. Laboratory results show that the 
formation efficien cy of fullerenes de pends on the content of hydrogen and this has led to 
the suggestion by lCami et al.l (120101 ) that under hydrogen-poor environment, the formation 
of Ceo is favored, and otherwise the chemical pathway favors the formation of Polycyclic 
Aromatic Hydrocarbon (PAH) molecules. This scenario is supported by the fact that the 
infrared spectrum of Tc 1 does not exhibit the aromatic infrared bands (AIBs) commonly 
assigned to as arising from PAH molecules. However, c ontrary to this scenario. Cm w as de- 
tected in hydrogen- containing PNs showing the AIB s (IGarcia-Hernandez et al.ll2010l ). The 
detection of C 60 in NGC 7023 (jSellgren et al.ll2010l ). an object with strong AIB emissions, 
also confirms that Cro and th e AIB carriers can coexist. As an alternative explanation, 



Garcfa-Hernandez et al.l (120101 ) suggested that both PAH molecules and fullerenes are formed 
from the destruction of hydrogenated amorphous carbon (HAC), and the non-detection of 
AIBs in Tc 1 is due to the longer survival time of Cgo molecules. A search for Cm in a large 



sampl e of R Coronae Borealis (R Cor Bor) stars was performed by I Garcia- Hernandez et al. 
(120 111 ). Ceo is clearly detected in one of the R Cor Bor stars showing AIBs, but Ceo is absent 



in the most hydrogen-poor ones. These results further support the premise that fullerene 
and the carrier of the AIBs can coexist. 

PPNs are the des cendants of the asymptotic giant branch (AGB) stars and the immedi- 
ate precursors of PNs (lKwokill993l ). The changes of infrared spectroscopi c properties in AGB 
stars, PPNs, and PNs have been known as a consequence of evolution (IKwokl 120041 ) . Thus 
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PPNs can provide a unique opportunity to investigate the formation history of Ceo in cir- 



cums tellar envelopes. We have obtained infrared spectra of a sample of PPNs (jZhang et al. 



20101 ) which cover the wavelen gths of four C fin vib rational transitions at 7.0 (z/ 25 ), 8.5 (z^), 
17.4 (z/ 27 ), and 18.9 /mi (f 28 ) fjNemes et al.lll994l ). These spectra therefore can serve as a 
platform to search for Ceo m PPNs. 

In this paper we report a new detection of C 60 in a PPN, IRAS 01005+7910. Unlike 
the other sources in our sample, IRAS 01005+7910 does not show the 21-/zm feature and its 
central star ha s a higher tempera t ure (~ 21, 5 00 K), suggesting that it a PPN about to enter 
the PN stage (jZhang et al.1120101 ). iHu I ( 120021 ) classified it as a B2 Ie star with V magnitude 
of 10.85. Its hydrog en B aimer lines show P Cygni profiles. Through a study of the high- 
resolution spectrum, iKlochkova et al.l (120021 ) concluded that it is a carbon-rich PPN with a 
luminosity log(L/L ) = 3.6 at a distance about 3kpc. 



Data 



The study makes use of the infrared spectra retrieved from the Spitzer Heritage Archive 
(SHA) 0. The observations were conducted with th e Infrared Spectrogr aph (IRS; Houck et 
al. 2004) on the Spitzer Space Telescope (Spitzer, IWerner et al.l 120041 ) in 20 04 and 2006. 
We ha ve carried out a systematic search for Ceo i n the ten PPNs studied by IZhang et al. 



( 120101 ) and found that among the studied sample IRAS 01005+7910 is the unique one that 
clearly exhibits the Ceo features. The spectra of IRAS 01005+7910 was obtained with the 
short-wavelength low-resolution module (SL; 5-14.5 /im) as part of the program 30036 (PI: 
G. Fazio), and the short- wavelength high-resolution module (SH; 9.5-19.5/im) as part of 
the program 9 3 (PI: D. Cruikshank) . Details of the dat a processing have been described 



elsewhere (e.g. iHrivnak et al.ll2009l ; ICerrigone et al.ll2009l ) and are not repeated here. How- 



ever, no observation utilizing long-wavelength high-resolution module (LH; 18.7-37.2 /im) 
was made for IRAS 01005+7910. 



3. Results 

Figure [1] shows that the Spitzer /IRS spectrum of IRAS 01005+7910 is dominated by a 
thermal dust continuum, the 11.5 /mi SiC emission, the aromatic infrared bands (AIBs) at 
6.2, 7.7-7.9, 8.6, 11.3, and 12.7/xm, and the 15-20 /mi plateau feature. The AIB and the 



1 http://sha. ipac.caltech.edu/applications/Spitzer/SHA/ 
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(2009 


) and 


Zhang et al. 


(2010) 



i this paper, we report the detection of the C 60 features 
at 7.04 ± 0.05, 17.4 ± 0.05, and 18.9 ± 0.04 /an. The fourth expected C 60 feature at 8.5 /im 
feature is badly blended with the AIB 8.6 /mi feature. These C 60 features have a width of 
0.31 ± 0.05 /im, much broader than the spectral resolution. The measured widths of the C 



features are similar to those seen in other PNs (see Table 1, iGarcfa- Hernandez et al.ll2010l ). 
These C 60 features are not found in any other PPNs in our sample. This suggests that these 
features likely share a common origin, and thus strengthens the identification of Ceo as their 
carrier. 

The continuum was fitted using the feature-free spectral regions, and was subtracted 
from the observed spectrum. As the features in PPN spectra are usually broad and blended 
with each ot her, we cond u cted a spectral decomposition using the IDL package PAHFIT 
developed by ISmith et al.l (120071 ) in order to accurately measure the feature fluxes. Drude 
profiles are assumed for the AIBs and Ceo features. The actual profile of the 15-20 /im 
plateau, where Ceo 17.4 and 18.9 /im features are superimposed, is poorly known. We have 
assumed two broad Gaussian profiles having a width of 1.3 /im and peaked at 16.1 and 
17.5 /im for the plateau (this causes only slight uncertainty in the flux measurements since 
the Cgo features are much narrower than the plateau). For the fitting, we have taken into 
account the AIBs at 6.2, 6.4, 6.6, 6.8, 7.4, 7.6, 7.9, 8.3, 8.6, and 16.5 /mi. Figures [2] 
and [3] give the zoomed-in view of these Ceo features and the fitting results. From the 
fitting results, we derived fluxes of (3.0 ± 0.3) x 10" 15 W/m 2 , (1.2 ± 0.5) x 10~ 15 W/m 2 , and 
(2.9 ±0.5) x 10 -15 W/m 2 for the Cqq 7.0, 17.4, and 18.9 /im features, respectively. The errors 
of the fluxes w e re est imated using the full covariance matrix of the least-squares parameters. 
Sellgren et al.l (120101 ) found that the Ceo 17.4 /im transition in the spectrum of NGC 7023 
is partially blended with an AIB feature. If this is also the case for IRAS 01005+7910, the 
estimated strength of the 17.4 /im feature would be the upper limit. 



Hrivnak et all (I2000f ) presented the ISO-SWS spectrum of IRAS 01005+7910 covering 
a wavelength range of 2.4-45.4/mi. Due to the lower sensitivity of the ISO SWS compared 
to the IRS, the C 6 o features are completely overwhelmed by the noise in the ISO spectrum. 
However, a strong feature at 30 /im is det ected by ISO, and its presence was subse quently 



confirmed by the Spitzer/IRS sp ectrum (ICerrigone et al 



seems to support the finding by iGarcfa-Hernandez et al 



20091; IZhang et al.l 120101 ) . This 



(120101 ) that all the Ceo sources 



exhibit the 30 /im feature. However, this correlation only applies to circumstellar sources 
as the Spitzer archive spectra of the reflection nebulae NGC 2023 and NGC 7023 (program 
40276, PI: K. Sellgren) do not show the 30 /im feature. The 30/im feature is comm only 
seen in carbon-rich AGB stars, PPNs, and PNs (IForrest et al 



has been attributed to solid magnesium sulfide (MgS; iGoebel fc Moseleyl Il985l ). However 



1981 



Volk et al.M 2002h. and 
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the identification of MgS as the carrier of the 30 fim feature is de batable as this feature is 
only detected in carbon-rich sources. Recently, IZhang et al.l (120091 ) found that the MgS dust 
mass in circumstellar envelopes is not enough to account for the observed feature strength. 
Therefore, carbonaceous compounds might be more likely to be the carrier of the 30 /im 
feature. In order to establish a connection between the 30 /im feature with the Ceo features, 
more Ceo sources need to be discovered. 

In Figure [TJ we co mpare the Spitzer/ 'IR S spectra of IRAS 01005+7 910 with two other 
C 60 sources, NGC 7023 JSelkren et alJboioh and Tc 1 JCami et alJboioh . All three sources 
have a strong infrared excess and IRAS 01005+7910 and Tc 1 have a very red (low color 
temperature) continuum. IRAS 01005+7910 and NGC 7023 show strong AIB features, which 
are absent in Tc 1. The fact that IRAS 01005+7910 does not show the narrow atomic lines 
as seen in the spectrum of Tc 1 is consistent with the object being a PPN and its central 
star is not hot enough to ionize the surround envelope. After subtracting the continuum, 
we found t h at th e spectral shape of IRAS 01005+7910 is similar to that of NGC 7023. 
Cami et al.l (120101 ) also detected a few weaker C70 features in Tc 1, which are not seen in 
IRAS 01005+7910 and NGC 7023. Assuming that all the sources have the same C 7 o/C 6 o 
strength ratio, the C70 features in IRAS 01005+7910 should be well below the detection 
limit. 



4. Discussion 



The relative intensities of the Cqq lines reflect the excitation conditions of the molecule. 



Assuming a thermal distribution of the vibrational states, ICami et al.l (120101 ) determined the 
excitation temperature of Tc 1 to be ~ 330 K, and suggested that the C ap molecules are in a 



solid. Similar temperature values we re also derived in th e PN s studied bvlGarcfa-Hernandez et al 
(120101 ). Using a similar procedure as lCami et al.l ( 120101 ) and iGarcia-Hernandez et al.l (120101 ). 
we have constructed a vibrational diagram for IRAS 01005+7910 from the observed fluxes 
of the three Cqq lines (Figure HJ). An excitation temperature of 460 ± 50 K is derived. 

The Ceo line ratios in IRAS 01005+7910 are I7.0//18.9 = 1-0 ± 0.3 and /17.4/I18.9 = 
0.4 ± 0.2, which are very cl ose to the values fou nd by in NGC 7023 (/7.0//18.9 = 0.82 ± 0.12 
and /17.4//18.9 = 0.42 + 0.02; Iselkren et alfaioh . Since the central star of IRAS 01005+7910 
has a similar temperature as NGC 7023 (~ 17, 000 K), it is possible t hat the C fi n mo l ecules 
in IRAS 01005+7910 and NGC 7023 are excited in a similar manner. ISelkren et all (boioh 
suggests that Ceo molecules in NGC 7023 are in the gas phase and they are excited by 
UV photons from the central star followed by radiative cascade. If this is also the case for 
IRAS 01005+7910, then the excitation temperature derived above is just a representation 
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of population distribution and is not related to other physical quantities such as kinetic 
temperature of the gas. 

Laborat ory measurements show that the wave lengths of gas-phase Ceo bands shift with 
temperature (IFrum et al. Hl99ll ; lNemes et al.lll994j ). The positions of the four infrared bands 
^25, ^26, ^27, and u 28 shifts from 6.97, 8.40, 17.41, and 18.82 at OK to 7.11, 8.55, 17.53, 18.97 
at 1000 K, respectively. Our measurements show that the wavelengths of all Ceo bands in 
IRAS 01005+7910 lie inside these ranges. Assuming that the frequencies of the C 60 bands 
have a linear dependence on the temperature, we can estimate the temperature as 100-600 K. 
The observed widths of the Ceo features are about 0.3 /im (see Section 3), which correspond 
to a wavenumber width of ~60 cm -1 and ~10 cm -1 for the 7 and 17.4/18.9 /im bands 
respectively. These values ca n be compared wit h the laboratory measured widths of about 



13 cm 1 of gas-phase C 60 (IFrum et al. Ill99ll ). suggesting that a gas-phase origin of the 



molecule cannot be ruled out. 



We estimate the abundance of Ceo following the same method of ISellgren et al.l (120101 ) 
by calculating the total strength ratios between the Cqq and AIB feature s and assuming 



that the carrier of the AIB features to contain 6 ± 2% of interstellar carbon (ICerrigone et al. 
20091 ). The strength ratio of Cqq to AIB emission is 0.01 in the observed wavelength range 
of IRAS 01005+7910, resulting in a percentage of ca rbon in Ceo of 0.06 + 0.02%. This value 



is slightl y lower than those obtained in NGC 7023 (ISelkren et aUBOlOh and the PN SMP 
S MC 16 (IGarcfa- Hernandez et al.ll2010l ). but a factor of 25 lower than that of Tc 1 estimated 
bv lCami et al.l (bold ). 



With the detection of Ceo m circumstellar envelopes, the next question is how they 
are formed. In the laboratory, Ceo can b e effectively produ ced from the vaporization of 
graphite in a hydrogen-poor environment. ICami et al.l (120101 ) proposed that fullerenes are 
prod uced only in hydrogen-p o or en velopes created by a late AGB thermal pulse. How- 



ever, iGarcfa-Hernandez et al.l (120101 ) argued that Ceo can be synthesized under hydrogen- 
containing environment. As shown in Figure. (U the spectrum of IRAS 01005+7910 ex- 
hibits both hydrogen- containin g; AIBs and C 6 o features, supporting the latter hypothesis. 
Garcfa-Hernandez et al.l (120101 ) also suggested that both C 6 o and PAHs are the products of 
decomposition of HACs. In this scenario, photochemical proce ssing can lead to dehydro- 
genation of the dust grains and form PAHs and Ceo molecules (IScott et al.lll997l ). On the 
other hand, the de hydrogenation of dus t grains can also induce the formation of H2 in the 
grain surfaces (e.g. iFleming et al. Il2010l ). The strongest H 2 line in the observed wavelength 
range is the -0 S(0) transition at 28 /xm. This line is no t detected in the sp ectrum of IRAS 



01005+7910 (IHrivnak et al.l l2000l : ICerrigone et all 12009 



have detected H 2 in two PPNs with no detectable Ceo (IZhang et al.ll2010l ). Because of this 



Zha ng et al.l 120101) . Moreover, we 



-7- 



lack of correlation between the presence of H2 and Ceo, we are unable to give addition support 
for the idea that the formation of Cgo is the result of dehydrogenation of HACs. 

Our detection suggests that fullerenes can be for med in the PPN stag e. So far, there is 
no definite detection of C 60 in AGB stars. Although IClayton et al.l (119951 ) noted a possible 
emission feature centered at 8.6 /^m in the spectrum of the bright AGB star IRC+10216 , 
the other Ceo features were not detected in the ISO-SWS spectrum (jCernicharo et al.lll999l ). 
The circu mstellar spectra o f AGB stars are dominated by silicates or silicon carbide emission 
features (IKwok et al.l 119971 1 . Although there are a small number of AGB stars exhibiting 
AIBs, the AIBs mainly emerge in the post-AGB phase. Is it possible t h at th e formation 
of Ceo is related to the emergence of the AIB features? ISellgren et al.l (120101 ) found the 
Cqq and AIB emissions in NGC 7023 have different spatial distributions and attributed this 
to the effect of UV-excitation. It can be argued that the Cqq and the AIB carriers are 
already present in the AGB phase of evolution but are not excited until the stars evolve to 
the PPN phase. However, comparisons between the spectra of AGB stars, PPNs, and PNs 
suggest a sequence of molecular synthesis, with acetylenes forming in the late AGB phase, 
leadi ng to the formation of diacetylenes, triacetylenes, and benzene in the PPN phase (IKwok 
20041 ) . Since these molecules are detected in absorption, the question of excitation does not 
arise. Since benzene is the first step toward the synthesis of aromatic materials, we can say 
with confidence that aromatic compounds only form after the AGB. If Ceo molecules are 
synthesized during the AGB, they should be detectable with absorption spectroscopy. 



5. Conclusions 



The detection of Cgo in a PPN as reported in this paper, together with the detection of 
this molecule in 5 PNs, confirm that the late stages of stellar evolution is a phase of active 
molecular synthesis. Beginning with simple diatomic molecules such as CO, CN, C2, dozens 
of gas-phase organic molecules have been seen in the stellar winds from AGB stars. The 
formation of acetylene during the late AGB phase is believed to lead to the formation of 
benzene in the post-AGB phase of evolution. This also coincides with the first detection 
of vibrational modes of aromatic and aliphatic compounds. From this study, we now learn 
that gas-phase Ceo molecules may also form during the same epoch. As the number of Ceo 
detection increases, we would be in better position to study the relationships between C 60 
and the carriers of other spectral features such as the AIBs, and 21- and 30-//m features. 



The detection of Ceo in the outflows from evolved stars also raises the possibility of 
the molecule being detected in the diffuse interstellar medium as the molecule is stable and 
should be able to survive journeys through the interstellar medium (jFoing fc Ehrenfreund 



-8 - 



1994 ). The fact that a large variety of presolar grains (e.g., SiC) have been detected in 



meteorites (jZinnerl Il998f ) raises the possibility that presolar Cgo can be incorporated into 



comets and asteroids and be detected in meteorites. In f act, C fi n and C 7 n, a s well as higher 



fullerenes, have been detected in the Allende meteorite (IBecker et al.lll999l ). The evidence 
for stellar synthesis of Ceo in the late stages of stellar evolution as presented in this paper 
therefore adds further support to the idea of chemical enrichment of the Solar System by 
stellar molecular products. 

We are grateful to Peter Bernath for helpful discussions. We also thank Kris Sellgren 
and Jan Cami for providing the IRS spectra of NGC 7023 and Tc 1 from the their papers 
(Sellgren et al. 2010, and Cami et al. 2010). This work is based on observations made with 
the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under a contract with NASA. This work was supported by a grant 
to SK from the Research Grants Council of the Hong Kong Special Administrative Region, 
China (Project No. HKU 7020/08P) and a grant to YZ from the Seed Funding Programme 
for Basic Research in HKU (200909159007). 
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Fig. 1. — The Spitzer/IRS spectrum of IRAS 01005+7910 compared to the IRS spectra of 
two other known Cqq sources (NGC 7023 from Sellgren et al. 2010, and Tc 1 from Cami 
et al. 2010). The vertical dotted lines mark the wavelengths of the Ceo lines at 7.0, 8.5, 
17.4, and 18.9 /xm. The dashed line represents a fit to the continuum. The other prominent 
features at 6.2. 7.7. 8.6. 11.3. a,nd 12.7 am are AIBs. 
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Fig. 2. — Continuum-subtracted Spitzer/IRS spectrum of IRAS 01005+7910 (black dots) 
between 6 and 8.8 /im. The expected positions of the C 60 features at 7.04 and 8.45 /zm are 
marked. The blue curves show the model fits to the individual AIBs (a list of the features is 
given in the text), with the total model spectrum given as the green curve. The red curve is 
a model fit to the 7.0 /im Ceo feature. The 8.5 //m Ceo line is badly blended with the 8.6 /im 
AIB feature and is not detected. 
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Fig. 3. — Continuum-subtracted Spitzer/IRS spectrum of IRAS 01005+7910 (black curve) 
between 15 and 19.5 /im. The positions of the Cgo lines at 17.4 and 18.9/im are marked. 
The blue curves show the model fits to the individual dust/plateau features, the red curves 
are the model fits to the 17.4 and 18.9 fim Cqo lines, and the total model spectrum is given 
by the green line. Note that the 15-20 /im plateau is fitted by two Gaussian profiles. 



- 14 - 




1000 



1500 
E u /k (K) 



2000 



Fig. 4. — Vibration diagram for the detected C 60 bands in IRAS 01005+7910. N u , g u , and E u 
are the population, statistical weight, and excitation energy of the upper level, respectively. 



